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Abstract 

A computational study of a distributed electric propulsion wing with a 40° flap deflection has been 
completed using FUN3D. Two lift-augmentation power conditions were compared with the power-off 
configuration on the high-lift wing (40° flap) at a 73 mph freestream flow and for a range of angles of attack 
from -5 degrees to 14 degrees. The computational study also included investigating the benefit of corotating 
versus counter-rotating propeller spin direction to powered-lift performance. The results indicate a large 
benefit in lift coefficient, over the entire range of angle of attack studied, by using corotating propellers that 
all spin counter to the wingtip vortex. For the landing condition, 73 mph, the unpowered 40° flap 
configuration achieved a maximum lift coefficient of 2.3. With high-lift blowing the maximum lift 
coefficient increased to 5.61. Therefore, the lift augmentation is a factor of 2.4. Taking advantage of the 
fullspan lift augmentation at similar performance means that a wing powered with the distributed electric 
propulsion system requires only 42 percent of the wing area of the unpowered wing. This technology will 
allow wings to be ‘cruise optimized’, meaning that they will be able to fly closer to maximum lift over drag 
conditions at the design cruise speed of the aircraft. 


Nomenclature 

Cy drag coefficient Vv freestream velocity, ft/sec 
CD HNL drag coefficient, high-lift nacelles contribution Vi ,atio ratio of tip speed to freestream velocity 
Cp pylons drag coefficient, pylons contribution WwW aircraft weight, lbf 
Corn drag coefficient, tip nacelles contribution (Xp, Vp» Zp) coordinates for the center of the propellers 
Cp,wing drag coefficient, wing contribution y axis along the wing span, in 
C, lift coefficient y nondimensional height, first node in boundary layer 
Cerf effective lift coefficient: C,+C, prop Symbols 
Ci max maximum lift coefficient a aircraft angle of attack, degrees 
Ch prop lift coefficient from the contribution of A delta 

propeller thrust in lift direction p density, slugs/ft° 
C sectional lift coefficient Abbreviations & Acronyms 
Ca pitching moment coefficient BSL Menter k-w basic turbulence model 
Cy pressure coefficient CFL Courant-Friedrichs-Lewy number 
Cre reference chord, ft. CW clockwise 
Cr thrust coefficient DEP distributed electric propulsion 
Cog torque coefficient HL high lift 
D drag force, lbf HLN high-lift nacelles, including pylons 
d propeller diameter, ft. KEAS knots equivalent airspeed 
h altitude, ft. LM Langtry-Menter transition model 
Kr normalized thrust coefficient max Maximum 
K, Q normalized torque coefficient mph miles per hour 
L lift force, Ibf prop propeller 
M Mach number QCR quadratic constitutive relation 
P pressure, Ibf/in? RPM revolutions per minute 
q dynamic pressure, Ibf/in SA Spalart-Almaras one equation turbulence model 
Re Reynolds number based on C,¢¢ SARC SA rotation and curvature correction 
S wing reference area, ft” hp horse power 
T temperature, °F SST Menter’s Shear Stress Transport model 
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I. Introduction 

NASA is taking the necessary first steps to reduce the environmental impact of general aviation 
aircraft, and potentially, larger aircraft as batteries continue to improve. In our future, NASA envisions 
electrically-powered airplanes filling the skies and eliminating the use of jet fuel that pollutes our 
atmosphere with carbon dioxide, nitrogen oxides, hydrocarbons and soot particles [1]. Researchers at NASA 
investigated distributed electric propulsion (DEP) in the Leading Edge Asynchronous Propulsion 
Technology (LEAPTech) project as a step to replace internal combustion engines with battery-powered 
motors that drive high-lift and cruise propellers. Aircraft with DEP systems will allow for a smaller wing 
(i.e., shorter span and higher aspect ratio), lower cruise drag, and consume 5 times less energy than a 
baseline, aviation gas-fueled airplane [2]. To assist takeoff and landing with a smaller wing, electrically- 
driven, high-lift propellers are mounted on nacelles upstream of the wing leading edge to increase dynamic 
pressure over the wing, increasing lift. The advanced-design propellers are designed to fold smoothly onto 
the nacelle during cruise flight. 

Since no previous DEP experiments had been conducted, a ground test experiment was designed with 
NASA SBIR funding. Figure 1 shows the DEP high-lift wing that was tested experimentally at the NASA 
Armstrong Flight Research Centeron the Hybrid-Electric Integrated Systems Testbed (HEIST) truck as part 
of the Leading Edge Asynchronous Propulsion Technology (LEAPTech) project. Initial computational fluid 
dynamics (CFD) results of the LEAPTech wing were published in reference 3. The current paper documents 
the results from a computational study using the FUN3D [4] Navier-Stokes CFD code for powered lift 
augmentation of the LEAPTech wing, assessing alternate approaches for the propeller spin rotation and 
various high-lift propeller combinations for blowing during powered lift. Although it would be desirable to 
compare the CFD results with the experimental data obtained from the HEIST test, the test data was not of 
sufficient quality to validate CFD data because of the uncertainty associated with this new truck-test 
technique. However, the LEAPTech test data was used to qualitatively compare with the CFD solutions to 
gain confidence that the tools could predict powered-lift increments comparable to those seen in the HEIST 
experiment. 


II. Configuration 

The configuration used in this report was the LEAPTech high-lift wing with a 40° flap deflection 
and DEP system that was tested on the HEIST truck. The wing consisted of the NASA GAW-1 airfoil with a 
fullspan 30% chord Fowler flap [5] and 18 propellers mounted on electric motors in nacelles upstream of the 
wing leading edge. The tip nacelles on the test article did not include the electric cruise motors that would be 
present to power the cruise propellers in a real airplane application. The LEAPTech wing has 4° twist and a 
= 0° is defined when the root section incidence is 5°. At the a = 0° inclination, all of the high-lift motor 
nacelles are at 0° incidence to the freestream flow. The wing reference area is 7944.48 in? (55.17 ft), 
reference span is 372 in., and the reference chord is 19.188 in. The moment reference center was located at 
the root quarter chord (6.9014 in., 0 in., 31.9185 in.). 


Ill. Computational Study 

Fully turbulent flow predictions of the LEAPTech wing were computed with the FUN3D Navier- 
Stokes CFD code using the Spalart-Allmaras (SA) one-equation turbulence model with rotation and 
curvature correction (RC) [6] and the quadratic constitutive relation (QCR) [7]. FUN3D is a node-based, 
finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. All of the FUN3D simulations 
for this study used an actuator disk representation for the high-lift propellers, with thrust and torque 
coefficient data used for inputs. The thrust and torque coefficient data were derived from the XROTOR [8] 
blade element momentum analyses. All cases were run with steady-state calculations that used a local time 
step pseudo-time advancement scheme. The pseudo-time advancement CFL number was ramped up from | 
to 100 over 1000 iterations for lower angles of attack and the maximum CFL was limited to 40 for the 
higher angles of attack. 

Two approaches for the propeller spin direction were investigated for blown-wing lift augmentation, 
to determine which approach generated the most lift. Figure 2(a) shows the counter-rotating propeller 
approach, which has pairs of propellers spinning toward each other. Figure 2(b) shows the propellers all 
spinning in the direction opposite the wingtip vortex in the corotating propeller approach. The propeller 
numbers are labeled in figure 2(b). Additionally, partially-blown wing studies were investigated to 
understand the impact of operating fewer propellers as well as the impact of inoperable propellers on lift. 
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The partially-blown configurations include corotating propellers 10 through 15 operating, corotating even 
propellers operating, and corotating propellers 11 through 16 operating. 


A. Grid Generation 

The GridTool/VGRID [9] software was used to generate the meshes for FUN3D simulations. 
GridTool is used to take a supplied IGES file and create the necessary points, curves, patches and grid 
topology (sources) to define the surface for grid generation. VGRID is a tetrahedral unstructured grid 
program that can either be run interactively or in a batch mode with grid spacing corresponding to the 
strength of user-defined sources placed in the domain. The advancing layers method [10] is used to generate 
the boundary layer and the advancing front is used to generate the farfield tetrahedral mesh. POSTGRID is 
used for additional post-processing, to close any open pockets that do not have cells and to improve the grid 
quality. 

An unstructured, mixed element mesh was used for the FUN3D simulations. The mixed element 
mesh had prisms in the boundary layer and tetrahedral cells in the farfield, with pyramids to transition 
between the boundary layer and farfield cells. The farfield boundaries were extended to approximately 100 
chord lengths away from the wing. A semispan computational model of the right half of the wing was used 
since all calculations were computed at 0° angle of sideslip. No special patching is required within the 
FUN3D grid for the propellers, except the clustering of cells to resolve the flow. The clustering is 
accomplished in VGRID with a cylinder volume source. 

A flat-plate, turbulent boundary layer calculator was used with reference chord length (c;¢,), 
Reynolds number based on C,¢¢, intended y’, and an intended number of node layers within the boundary 
layer, to determine the inputs for VGRID of the first node height and the expansion rate of the grid. As 
recommended in reference 11, the intended y* used was 0.67 for a medium mesh resolution, and the number 
of node layers in the boundary layer was specified as 26. The corresponding VGRID spacing for the first 
node height was 2.8e-4 in. The boundary layer specifications were good for the unblown wing because the 
solutions had y" of 1 or less. The y” values are not shown in this paper. However, for the blown-wing cases, 
the y” varied with values less than 1 in the unblown region and values less than two on the wing upper 
surface in the blown region. Along the wing leading edge, the y’ values were less than three on the ‘blade 
up’ side of the nacelle where the upwash increases the local velocity, and y* values were less than two on the 
‘blade down’ side of the nacelle. 

An initial grid study was completed with coarse, medium and fine meshes to determine an 
appropriate mesh for predicting C,. The coarse, medium and fine meshes had 35.4, 42.9, and 53.6 million 
total mesh points, respectively. The viscous boundary layers had 24.3, 29.3 and 35.9 million mesh points for 
the coarse, medium and fine meshes, respectively. The lift, drag and pitching moment coefficients for these 
meshes are shown in figure 3. The coefficients are similar for the fine and medium meshes, while the coarse 
mesh data is the outlier. The coarse mesh does not have enough grid resolution to adequately model the 
flow, predicting a severe, premature stall break that isn’t present in the medium and fine mesh data. The 
wing stalls after a = 2° for the coarse mesh with a severe reduction of lift of AC, = -0.34 (10.7%), an 
increase in drag of ACp= 0.14, and a AC,,= 0.35 reduction in negative pitching moment from a = 2° toa = 
4°. Both the medium and fine meshes predict maximum lift coefficient at a = 4° for the angles of attack 
studied, however, a finer a scale may refine C; maq,. The medium mesh grid was used for the remainder of 
this work to take advantage of fewer resources than required for the fine mesh grid. The y’ values were 
similar for all three meshes because the first node height did not vary between the meshes. In future blown- 
wing CFD, care must be used when choosing a first node height to account for the increased velocity from 
the propellers, above the freestream value, but not restrict the first node height to a value that is so small that 
it hinders convergence. 


B. Boundary Conditions and Propeller Modeling 

A symmetry plane boundary condition is used for the semispan grids. A no-slip boundary condition 
was used on all solid surfaces. Riemann invariants boundary condition was used for the farfield mesh 
boundaries. 

An actuator disk was used to simulate the high-lift (HL) propellers. There are 6 options to loading 
the rotor in FUN3D but the one implemented for this work was a body force based on the optimal 
distribution of Goldstein [12] implemented as described by Stern, Kim, and Patel [13]. Inputs for the 
actuator boundary condition include tip radius (8.79 in.), nacelle hub radius (2.88 in.), the location of the 
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nacelle hub centers for each propeller (x,, y,, Zp), the orientation of the propeller within the grid and the 
direction of rotation. Additionally, the ratio of tip speed to freestream velocity (Vi;atio), normalized thrust 
coefficient (K;) and normalized torque coefficient (Kg) were required for each propulsor and the values are 
shown in Table 1. The Vi ,q¢;9 was calculated with equation 1. The torque and thrust coefficients were 
calculated with equations 2 and 3 using the propeller diameter (d) of 1.465 ft. The normalized torque 
coefficient and thrust coefficient were computed with equations 4 and 5 and used in the actuator boundary 
condition input file. 


Table 1. Values for the actuator disk model to simulate HL propellers, h = 2300 ft. 


Description hp hp/prop | RPM V, mph T, °F K Q Kr V.,atio 
Windmilling 23.4 1.3 3000 | 73 40 0.03330 0.01541 2.1485 
Takeoff 257.4 14.3 6370 | 73 40 0.04555 0.01764 4.5621 
Grid and Propeller | 300.6 16.7 6147 | 73 60 0.05472 0.02377 4.3876 
Rotation Study 
Partial-Span varies 16.7 6,771 | 60 60 0.05846 0.02419 5.3251 
Blown Wing Study 
Landing 324 18 6860 | 73 40 0.04628 0.01778 4.9130 
RPM 
V, ratio = [n+ ( 60 ) «al /V (1) 
RPM 
Co =Torque/[pC7a]_— 2) 
Cr = Thrust /[pC—)7d*| (3) 
¢ 4 
Kg =3 (4) 
4 
Kp = 73 Cr (5) 


C. Flow Conditions 

All of the solutions were computed at an altitude of h = 2300 ft. to best match conditions at the 
NASA Armstrong Flight Research Center. 

The initial computational grid study and propeller rotation study were computed with a freestream 
velocity of 73 mph, M = 0.096, an air temperature of T = 60°F, and Re = 999,569 (1 million). The partially- 
blown wing study was conducted with a freestream velocity of 60 mph, M = 0.079, an air temperature of T = 
60°F, and Re = 818,503 (0.82 million) for various combinations of corotating propeller blowing for lift 
augmentation. The power setting for the grid study and the propeller rotation study was 300.6 hp (16.7 
hp/prop) at 6,147 RPM. The power setting for the partially blown wing study was 16.7 hp/prop at 6,771 
RPM. 

The data in figures 10 through 13 had freestream conditions of 73 mph, M = 0.098, angles of attack 
from -5° to 14°, T = 40°F, and Re = 1,069,038 (1.07 million). The 73 mph speed at NASA Armstrong, 
corresponds to a 61 knots equivalent airspeed (KEAS) stall speed at sea level and also was close to the 
maximum possible speed of the HEIST truck. The solutions in figure 10 were conducted to compare with 
the test propeller RPM and ambient conditions of the HEIST test at AFRC. Two powered RPM settings 
were used that would correlate with takeoff and landing power requirements (Table 1). The high-lift 
propellers were permanently fixed upstream of the wing and were not folding, so it was not possible to 
experimentally test a power-off case. A third RPM setting was tested that represented a windmilling 
condition (Table 1), with just enough input power so the propellers were being driven forward by the motors 
and not driven by the oncoming flow. 
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D. Convergence Criteria 

The criteria used to monitor and determine solution convergence was a drop of at least two orders 
of magnitude for the flow solution residual. In addition, for an interval of 2,000 iterations, a solution was 
finished running if the standard deviation of C,; and Cp were less than 0.0020. The unblown cases all 
converged well, with a standard deviation of less than 0.00014 for C,, 0.00002 for Cp, and 0.00015 for C,,. 
The windmill-power blown cases all converged well within the stated criteria. The takeoff-power blown 
cases all converged within the stated criteria for a < 4° but the standard deviation for C, ranged from 0.0026 
to 0.0071, from 0.0015 to 0.0034 for Cp, and from 0.0026 to 0.0105 for C,, at a = 4°. The landing-power 
blown cases all converged within the stated criteria for a < 8°, but the standard deviation for C, ranged from 
0.0026 to 0.0046, from 0.0020 to 0.0033 for Cp, and from 0.0053 to 0.0091 for C,,, at a=10° and a =11°. 
These standard deviations at the high angles of attack with power were slightly higher than planned but were 
deemed acceptable near wing stall. 


IV. Results 

Figure 4 shows comparisons of lift, drag, and pitching moment coefficients for the two propeller 
spin approaches, corotating and counter-rotating, for the blown, high-lift wing (40° flap). The conditions are 
73 mph and a power setting of 300.6 hp (16.7 hp/prop) at 6147 RPM. Although the HEIST truck test was 
originally designed with counter-rotating propellers, these computational results indicated a larger benefit in 
lift coefficient over the entire range of angle of attack by using corotating propellers (fig. 4(a)). The 
maximum lift coefficient occurs at a = 4° for both rotation approaches, but the corotating propellers increase 
lift coefficient by AC, = 0.31 compared to counter-rotating propellers. The drag coefficient is also increased 
by a maximum of ACp = 0.11 over the range of angle of attack (fig. 4(b)), but the high-lift propellers are 
only used during take-off and landing operations. The advanced-design propellers will be designed to fold 
smoothly onto the nacelle during cruise flight. The pitching moment coefficient is more negative for the 
corotating propeller approach as compared to counter-rotating propellers, but the trends of C,, with angle of 
attack are similar (fig. 4(c)). 

Spanloading distributions of section lift coefficient (c,) along the span (y) can be used to further 
evaluate the lift coefficient data in figure 4(a). The effect of angle of attack on spanloading for the counter- 
rotating propeller approach and the corotating propeller approach are shown in figures 5(a) and 5(b), with 
nacelle numbers above the spanwise location of the nacelle hubs and directional arrows indicating the 
propeller rotation. The spanloading for the unblown (no power) wing is also shown in both figures to 
emphasize the benefit of blowing on lift coefficient, with unblown wing c,; curves well below the blown 
wing c; curves from y = 22 to y = 182, between nacelles 10 to 18. The blown wing c, distributions for the 
counter-rotating propeller approach (fig. 5(a)) for a = 0° and a = 4° are similar, but c, for a = 4° is slightly 
higher than a = 0° from 0 < y < 20 and at the hub of nacelles 15 and 16. The lift coefficient drops 
substantially at a = 8° and a = 9° as the flow separates behind nacelles 15 and 16, from the left outer tip of 
propeller 15 to the hub of nacelle 16, while the rest of the c; distributions are mostly similar. The blown 
wing c, distributions for the corotating propeller approach (fig. 5(b)) for a = 0° and a = 4° are also similar, 
but c, for a = 4° is slightly higher than a = 0° from 0 < y < 15 and at nacelle hubs 11 to 18. The lift 
coefficient drops at a = 8° and a = 9°, but not as much as the counter-rotating propeller approach, as the 
flow separates from y = 140 to y = 156, between nacelles 16 and 17. 

The effect of the propeller spin approach on spanloading is shown in figure 6 for angles of attack at 
0°, 4°, 8° and 9°. Propeller ‘blade up’ locations are defined as regions where the blade is moving upward at 
the wing leading edge, which is to the left side (pilot’s view) of nacelle 10 and between nacelles 11 and 12, 
13 and 14, 15 and 16, and 17 and 18 for the counter-rotating propeller approach, and to the left side (pilot’s 
view) of all the nacelles for the corotating propeller approach. Propeller “blade down’ locations are defined 
as regions where the blade is moving downward at the wing leading edge, which is between nacelles 10 and 
11, 12 and 13, 14 and 15, 16 and 17 and to the right side (pilot’s view) of nacelle 18 for the counter-rotating 
propeller approach, and to the right side of all the nacelles for the corotating propeller approach. In figure 6 
for angles of attack 0° and 4°, the counter-rotating propeller approach has higher c, in the ‘blade up’ regions 
due to upwash from two propellers, rather than from one propeller in the corotoating approach. However, the 
counter-rotating propeller approach has much lower c, than corotating propellers in the ‘blade down’ regions. 
Again, at angles of attack of 8° and 9°, the counter-rotating propeller approach has a separated flow region 
between approximately y = 106 in. to y = 140 in., which is between the inboard tip of propeller 15 and 
nacelle hub 16 (figs. 6(c) and 6(d)), which results in a large loss of lift (fig. 4(a)), an increase in drag (fig. 
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4(b)), and a shift in pitching moment (fig. 4(c)). The separated flow region for the corotating propeller 
approach is slightly more outboard, at a span location between nacelle hubs 16 and 17, than counter-rotating 
propeller approach, and the c, in this region is still as high as the counter-rotating propellers (figs. 6(c) and 
6(d)), so the loss of lift is not as severe for the corotoating propeller approach. 

To further visualize the effects of blowing from each propeller spin approach, the pressure 
coefficient on a plane through the wing three inches back from the leading edge, and on the upper surface of 
the wing, for the maximum lift coefficient case at a = 4°, is shown in figures 7 and 8, respectively. The 
white regions in figure 7 represent the wing and nacelles geometry. The nacelles are labeled 10 to 18 in both 
figures. The regions of “blade up’ are easily visualized with red contours between pairs of nacelles for the 
counter-rotating propeller approach in both figures. The increased lift with the corotating propeller approach 
results from higher negative pressure coefficient values over the full span of the wing upper surface (fig. 
8(b)), as opposed to regions of higher negative pressure coefficient only occurring between the pairs of 
counter-rotating propellers (fig. 8(a)). There is a coupling of alternating upwash and downwash changes 
between pairs of propellers with the counter-rotating propeller approach. Based on the improved lift with 
corotating propellers, the project decided to proceed with the corotating propeller approach for the remainder 
of the HEIST truck test, and with future wing designs. 

A partial-span blowing scenario study was conducted. The LEAPTech wing had full-span blowing 
implemented. It was known that the X-57 Maxwell configuration would not have blowing in the wing-tip 
region, so the performance of inboard blowing only was of interest. As a failsafe mechanism, blowing 
motors can be installed on a double bus electrical system. It was of interest to document the degradation in 
high-lift augmentation if one electrical bus would fail. Finally, the sensitivity to the spanwise location of the 
blowing region was investigated. Figure 9 shows the comparisons of upper-surface pressure coefficient for 
the partially-blown wing study at 60 mph, a = 4°, and a power setting of 16.7 hp/prop at 6,771 RPM for 
each of the corotating propellers that are operating. In this study, all of the powered-lift configurations were 
propellers corotating in a direction opposite of the wing-tip wake vortex. A table of lift, drag, and pitching 
moment coefficients is also shown for the following blown configurations: (a) all propellers (10-18) 
operating, (b) partially-blown wing with propellers 10 through 15 operating, (c) partially-blown wing with 
even-numbered propellers operating, and (d) partially-blown wing with propellers 11 through 16 operating. 
The case where all of the propellers are operating generates the most lift with C, = 6.422 and pressure 
coefficient reaching C,, of -13 along the full span in the blown-wing region (fig. 9(a)). The cases where three 
of the propellers are not operating, 10 to 15 corotating (fig. 9(b)) or 11 to 16 corotating (fig. 9(d)), both 
achieve a HL coefficient of C, = 5, but reduced from C, = 6.422 as the blown-wing region shrinks. The drop 
in C, compared to all propellers operating is a result of less blowing along the span and a smaller region of 
pressure coefficient values reaching at least C;, = -13 along the leading edge. The even-numbered propellers 
operating still have a HL coefficient of C, = 4.574 with regions of upwash on the inboard side of the even- 
numbered propellers creating a C,, peak suction of -13 (fig. 9(c)), but lift is substantially lower than 
operating all the propellers. The partially blown wing results in 22 to 29 percent reduction in lift compared 
to the fully-blown wing with all propellers operating. These results emphasize the need for system 
redundancy to avoid a substantial penalty in lift from a partially blown wing due to motor or propeller 
failure. In the partially-blown wing case with even-numbered propellers operating, the lift is reduced by 
29%. This means that compared to the fully-blown wing case, if the original stall speed was 61 KEAS, then 
a partially-blown wing case with even-numbered propellers operating would only have a stall speed of 72 
KEAS. Future studies might investigate the additional power required for the partially-blown wing operable 
propellers to achieve the lift of the fully-blown wing with all propellers operable. 

Runs were conducted to compare with the test propeller RPM and ambient conditions of the HEIST 
test at the NASA Armstrong Flight Research Center. Two powered RPM settings were used that would 
correlate with a takeoff and landing power requirement. The high-lift propellers were permanently fixed 
upstream of the wing leading edge on the nacelles and were not folding, so it was not possible to 
experimentally test a power-off case. A third RPM setting was tested that represented a windmilling 
condition, with just enough input power so the propellers were being driven forward by the motors and not 
driven by the oncoming flow. A comparison of lift, drag and pitching moment coefficients between unblown, 
windmill power, takeoff power, and landing power for the high-lift wing with corotating propellers at 73 
mph is shown in figure 10. The maximum lift coefficient unblown is C; mg, = 2.318 at a = 1°. The 
maximum lift coefficients of Cy may = 5.252 at takeoff power (257.4 hp) and C; max = 5.611 at landing 
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power (324 hp) both occurred at a = 4°. Therefore, the maximum lift augmentation, C, 4, from the 
unblown wing (2.318) to the blown wing (5.611), is a factor of 2.42. The blowing lift augmentation does 
come with an increase in drag (fig. 10(b)), from Cp = 0.3209 to Cp = 1.0443 at a = 4°, for example. The 
blowing lift augmentation also comes with a large impact on pitching moment (fig. 10(c)), from C,, = -1.720 
to Cy, = -4.894 at a = 4°, for example. In general, there is a fairly small change in pitching moment 
coefficient at each power condition over the range of angle of attack for the data shown. Note the pitching 
moment reference center is located at the root chord and not the mean aerodynamic chord of the wing. 

A comparison of the pressure coefficient contours between unblown, blown with takeoff power, 
and blown with landing power is shown in figure 11 at 73 mph and a = 4°. The contours shown on a 
spanwise plane through the wing is located 8 inches back from the wing leading edge. This figure shows 
much higher negative pressures over the wing for blown conditions that cause increased dynamic pressure, 
which is the benefit of the DEP system. Without lift augmentation (unblown), the lowest pressure coefficient 
above the wing in the spanwise plane is approximately C, = -1.5 (fig. 11(a)). With lift augmentation, the 
region of suction above the wing in the spanwise plane is substantially increased over most of the wingspan 
(out to nacelle 16) at takeoff power (fig. 11(b)), and to a further spanwise extent (out to nacelle 18) for the 
landing power (fig. 11(c)). Figure 12 shows the comparison of surface pressure coefficient between the 
unblown wing and the two blown-wing conditions at various span locations along the high-lift wing at 73 
mph and a = 4°. Figure 12 also shows the locations of the span stations on the wing (with black text) and the 
nacelles are numbered in red. The pressure coefficient plots further show the lift augmentation illustrated by 
the area between the higher negative pressures for the blown wing (red and blue curves) and the lower 
negative pressures on the unblown wing (green curve). It is interesting that at eight of the ten stations shown, 
the peak pressures on the leading edge are greater for the takeoff power (257.4 hp) than the landing power 
(324 hp), because the opposite is expected due to a greater lift coefficient for the landing power condition. 
However, there is more lift generated over the aft upper wing surface and on the flap for the landing power 
(red), compared to the takeoff power (blue), as shown by higher negative pressures. One explanation could 
be that the flow through the propeller gets straightened, effectively lowering the angle of attack, but the 
higher power setting yields a higher dynamic pressure (q) above the wing to generate more lift. Figure 12(h) 
shows the difference in pressures just outside of nacelle 16 at y = 145.8 in. where the takeoff power solution 
has a much lower and flat C;, distribution than the landing power that was also viewed in figures 11(b) and 
11(c). The unpowered condition looks separated downstream of the leading edge at y = 145.8 in. (fig. 12(h)) 
and there appears to be an overall loss of lift in the spanwise region between nacelles 16 and 17 in the 
pressure contour plot (fig. 11(a)). 

Taking advantage of the full-span lift augmentation at similar takeoff and landing performance 
means that a wing powered with the DEP system could be 42% of the wing area of the unpowered wing. 
This technology will allow wings to be ‘cruise optimized’, meaning that they will be able to fly closer to the 
maximum lift to drag ratio (L/Dingx) conditions at the design cruise speed of the aircraft. A comparison of 
spanwise lift coefficient between maximum landing power (18 hp/prop) and no power, for the high-lift wing 
(40° flap) is shown in figure 13. The benefit of DEP is clear in the increased sectional lift coefficient with 
high-lift blowing power, as compared to the lift coefficient of the unblown wing. Sectional lift coefficients 
greater than 7.5 are computed in the blowing region. The peaks and valleys in the powered sectional lift 
coefficient curve results from the propellers creating areas of increased upwash and downwash across the 
span that result in higher and lower lift on the wing. 


V: Conclusions 

A computational study of a distributed electric propulsion high-lift wing with a 40° flap deflection 
has been completed using FUN3D. The distributed electric propulsion wing in this study was tested 
experimentally on the Hybrid-Electric Integrated Systems Testbed truck at the NASA Armstrong Flight 
Research Center. 

The computational study included three lift-augmentation power conditions of the high-lift wing in 
73 mph freestream flow, in addition to the unblown wing without lift augmentation. The angle of attack 
ranged from -5 degrees to 14 degrees for various portions of the computational study. 

The computational study also included investigating the spin approach of the propellers, corotating or 
counter-rotating. The corotating propeller approach has all of the propellers spinning in a direction that 
counters the wingtip vortex, while the counter-rotating propeller approach has pairs of propellers spinning 
toward each other. The computational results indicate a larger benefit in lift coefficient over the entire range 
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of angle of attack, from -6 degrees to 9 degrees, by using corotating propellers. A maximum lift coefficient 
occurred at 4 degrees angle of attack for both rotation approaches, but the corotating propeller approach 
increased lift coefficient by 0.31 at 4 degrees angle of attack. 

At 73 mph, the high-lift wing with 40° flap deflection achieves a maximum lift coefficient of 2.3 
without high-lift blowing. Using the corotating propeller approach, the maximum lift coefficient is 5.25 with 
high-lift blowing at takeoff power, and 5.61 with high-lift blowing at landing power. Therefore, the lift 
augmentation is a factor of 2.4. The blown-wing lift augmentation does come with an increase in drag, 
however, the lift augmentation occurs during takeoff and landing, and thus, drag is improved during cruise 
when the high-lift props are not operating and fold smoothly onto the nacelles as intended. The blowing lift 
augmentation also comes with a large impact on pitching moment but there is a fairly small change in 
pitching moment coefficient for each power condition over the range of angle of attack. 

Taking advantage of the fullspan lift augmentation at similar performance means that a wing 
powered with the DEP system could be 42% of the wing area of the unpowered wing. This technology will 
allow wings to be ‘cruise optimized’, meaning that they will be able to fly closer to ZL/Ding. conditions at the 
design cruise speed of the aircraft. 
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f _— \ s — e 
Figure 1. The DEP high-lift wing tested at NASA Armstrong on the Hybrid-Electric Integrated 
Systems Testbed (HEIST) truck. 


(a) Counter-rotating propellers on the right semispan wing 
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(b) Corotating propellers on the right semispan wing 
Figure 2. Illustration of two approaches for the propeller spin direction. 
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(a) Lift coefficient 
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(c) Pitching moment coefficient 


Figure 3. Grid study for the high-lift wing (40° flap) at 73 mph, M = 0.096, Re = 1 million, h = 2300 ft, 
T = 60°F, and 300.6 hp (16.7 hp/prop, 6147 RPM). 
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(c) Pitching moment coefficient 
Figure 4. The comparison of lift, drag and pitching moment coefficient between propeller spin 
direction approaches for the blown, high-lift wing (40° flap) at 73 mph, M = 0.096, Re = 1 million, 
h = 2300 ft, T = 60°F, and 300.6 hp (16.7 hp/prop, 6147 RPM). 
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(a) counter-rotating propeller approach 
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(b) corotating propeller approach 
Figure 5. The effect of angle of attack on spanloading for each propeller spin direction approach, for 


the blown, high-lift wing (40° flap) at 73 mph, M = 0.096, Re = 1 million, h = 2300 ft, T = 60°F, and 
300.6 hp (16.7 hp/prop, 6147 RPM). 
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Figure 6. The effect of propeller spin direction approach on spanloading at four angles of attack for 
the blown, high-lift wing (40° flap) at 73 mph, M = 0.096, Re = 1 million, h = 2300 ft, T = 60°F, and 
300.6 hp (16.7 hp/prop, 6147 RPM). 
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Figure 6 Continued. The effect of propeller spin direction approach on spanloading at four angles of 
attack for the blown, high-lift wing (40° flap) at 73 mph, M = 0.096, Re = 1 million, h = 2300 ft, T = 60°F, 
and 300.6 hp (16.7 hp/prop, 6147 RPM). 
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(a) Counter-rotating propellers 


(b) Corotating propellers 
Figure 7. The comparison of pressure coefficient on a spanwise plane between approaches to propeller 
spin direction for the blown, high-lift wing (40° flap) at 73 mph, M = 0.096, a = 4°, Re = 1 million, h = 
2300 ft, T = 60°F, and 300.6 hp (16.7 hp/prop, 6147 RPM). 
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(a) Counter-rotating propellers 
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(b) Corotating propellers 
Figure 8. The comparison of pressure coefficient on the wing upper surface between approaches to 
propeller spin direction for the blown, high-lift wing (40° flap) at 73 mph, M = 0.096, a = 4°, Re =1 
million, h = 2300 ft, T = 60°F, and 300.6 hp (16.7 hp/prop, 6147 RPM). 
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Figure 9. The pressure coefficient on the upper wing surface at 60 mph, M = 0.079, a = 4°, Re = 0.82 
million, h = 2300 ft, and T = 60°F for various corotating propellers blowing for lift augmentation using 
16.7 hp/prop at 6,771 RPM. 
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(c) pitching moment coefficient 
Figure 10. The lift, drag and pitching moment coefficients for the high-lift wing (40° flap) at 73 mph, 
M = 0.098, Re = 1.07 million, h = 2300 ft, and T = 40°F. The blown wing was computed with 3 power 
settings using corotating propellers. 
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(c) Corotating, landing power 324 hp (18 hp/prop, 6860 RPM) 
Figure 11. Comparison of pressure coefficient between the unblown wing, the blown-wing with 


corotating takeoff power, and the blown-wing with corotating landing power high-lift wing (40° flap) 
at 73 mph, M = 0.098, a = 4°, Re = 1.07 million, h = 2300 ft, and T = 40°F. 
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Figure 12. Pressure coefficient for the high-lift wing (40° flap) at 73 mph, VM = 0.098, a = 4°, Re = 1.07 


million, h = 2300 ft, and T = 40°F. 
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Figure 12 Continued. Pressure coefficient for the high-lift wing (40° flap) at 73 mph, M = 0.098, a = 4°, 
Re = 1.07 million, h = 2300 ft, and T = 40°F. 
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Figure 13. Comparison of spanwise lift coefficient between the unblown high-lift wing (40° flap) and 
the blown-high-lift wing (40° flap) with corotating propellers at landing power (324 hp, 18 hp/prop, 
6860 RPM) at 73 mph, M = 0.098, a = 4°, Re = 1.07 million, h = 2300 ft, and T = 40°F. 
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